Decoherence and mode-hopping in a magnetic tunnel junction-based spin-torque 
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We discuss the coherence of magnetic oscillations in a magnetic tunnel junction-based spin-torque 
oscillator as a function of external field angle. Time-frequency analysis shows mode-hopping between 
distinct oscillator modes, which arises from linear and nonlinear couplings in the Landau-Lifshitz- 
Gilbert equation, analogous to mode-hopping observed in semiconductor ring lasers. These couplings 
and therefore mode-hopping are minimized near the current threshold for antiparallel (AP) align- 
ment of free layer with reference layer magnetization. Away from the AP alignment, mode-hopping 
limits oscillator coherence. 
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Magnetization precession at GHz frequencies can 
be sustained in spin-valve (SV) and magnetic tunnel- 
junction (MTJ) based spin-torque oscillators (STOs) by 
directly transferring spin angular momentum from spin 
polarized current [1, 2 to the free-layer magnetization or- 
der parameter. A fundamental question is the stochastic 
phenomena that govern its coherence time (r c ). Theoret- 
ical studies [3j [4] have investigated decoherence through 
thermal noise, assuming that only a single mode is ex- 
cited. Other theoretical works conclude that only the 
lowest energy mode supports sustained oscillations [5]. 
Yet experiments clearly show the existence of multiple 
modes in SVs [6H9] and MTJs [T0HT4] , and persistent 
mode-hopping [8j [9] between several modes. The impact 
on r c of such mode-hopping is largely unexplored, and 
theoretical investigations of its origin are entirely lack- 
ing. 

Here we present systematic experimental investigations 
into mode- hopping and its impact on coherence time, 
as functions of current and applied field angle in MTJ- 
STOs. We derive equations showing how such systems 
are analogous to semiconductor ring lasers (SRLs), and 
under driving can exhibit mode-hopping in the pres- 
ence of stochastic noise. Nonconservative torques in 
the Landau-Lifshitz equation couple individual modes. 
The coupling has a minimum when the nonconservative 
torques cancel each other, explaining the experimentally 
observed suppression of mode-hopping at angles near 
(but not exactly at) antiparallel (AP) alignment. Mi- 
cromagnetic simulations of magnetization dynamics sup- 
port this picture. Finally, we show that although mode- 
hopping is the limiting factor for STO coherence at most 
angles, single-mode nonlinear spin-torque auto-oscillator 
(NSTO) theory [3] qualitatively holds when dwell-time is 
sufficiently large compared with the oscillations' coher- 
ence time. 

The MTJ nanopillars used in this work are similar to 
those in Ref. 14. The layer structure consists of IrMn 
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Figure 1. (color online). (a) Experiment (circles) and calcu- 
lated (solid line) resistance versus in-plane field angle ip at 
H=450 Oe; ip — 180° corresponds to AP alignment of FL 
and RL. Inset: magnetoresistance loop measured at ip = 180°. 
(b) Map of power (dB) vs. frequency (/) and ip for 7=8 mA 
and #=450 Oe. 



(5)/CoFe (2.1)/Ru (0.81)/CoFe (l)/CoFeB (1.5)/MgO 
(l)/CoFeB (3.5) (thickness in nm), where the bottom 
CoFe layer is the pinned layer (PL), the composite 
CoFe/CoFeB represents the RL, and the top CoFeB layer 
is the FL. We discuss results from a circular device with 
approximate diameter 240 nm, resistance-area product 
1.5 Q /im 2 , and tunneling magnetoresistance 75%. The 
RL magnetization equilibrium direction is along the pos- 
itive x-direction, which is also 0° of the applied field. We 
use the convention that a positive current flows from the 
FL to the RL. 

Figure Ilia) shows the resistance R as a function of 
the in-plane angle ip of applied field H. Here we fo- 
cus on magnetic excitations with H = 450 Oe and 
140° < <p < 220°, for which the FL rotates coherently 
with the field [Fig. Ilia)]; with a positive current, FL 
modes are excited. In general, several modes are found 
in frequency-domain measurements, and we have identi- 
fied five potential FL modes (m^, i = 1, . . . , 5) [14 from 
measurements at 8 mA. These modes' frequencies de- 
crease when the field angle moves away from 190°, and 
are asymmetric about their maxima. Both the decrease 
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Figure 2. (color online). Wigner distribution (first column), 
and normalized autocorrelation envelope (second column) at 
(a, b) 140°, (c, d) 196°, and (e, f) 220° for I = 8 mA and 
H — 450 Oe. Red lines are exponential fits to experimental 
autocorrelation envelope (dashed lines). 



and asymmetry are unexpected from a simple Kittel for- 
mula [14 J, in which ferromagnetic coupling between the 
FL and RL would result in mode frequencies symmetric 
about 180° and frequency minima at 180°. 

In order to quantitatively analyze the STO's time- 
frequency behavior, we collected 1,000 ns-long time- 
traces [15] of the STO voltage signal at different field 
angles and bias currents. Figure [2] shows the frequency 
vs. time obtained from Wigner transforms of the traces 
at three representative angles, cp = 140°, 196°, and 220°, 
and at 8 mA using a time (frequency) resolution of 1.25 ns 
(0.8 GHz). We chose these angles because at 140° and 
220° the FL and RL are far from their AP configuration, 
but the FL still rotates coherently with the external field, 
and at 196° we observed longest coherence time of mode 
ni2, because of the device's slight asymmetry [Fig. [lib)] 
and because closer to 180°, the observed ni2 coherence 
time is limited by frequency doubling [TTJ [16] . The plots 
in Fig. [2] show that for 196°, the STO mostly stays in 
mode ni2. In contrast, frequent mode- hopping occurs be- 
tween mi and m2 at 220°, and between mi, ni2, and 1114 
at 140°. We note that this mode-hopping behavior be- 
tween distinct STO modes differs from the observations 
in Ref. [17] where the dominant mode's frequency fluctu- 
ated in time, on the scale ~ 0.01 GHz. Detailed exami- 
nation of time-frequency plots shows that mode-hopping 
occurs at all field angles and currents [15], is completely 
random in time. The central mode ni2 is the most stable 
at all angles; spin torque preferentially excites mode m 2 . 

The coherence time of mode ni2 can be obtained by 
autocorrelating time-traces. The right of Fig. [2] shows 
the traces' normalized autocorrelation function, T(r), fil- 
tered in a range of -300 MHz/+400 MHz around m 2 . This 
filter width avoids overlapping other modes. The auto- 
correlation functions decay exponentially in time, consis- 
tent with thermally activated stochastic processes lead- 
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Figure 3. (color online). Coherence time r c (filled symbols) 
and average dwell-time t ave (open symbols) vs. I at 140° 
(black squares), 196° (red circles), and 220° (blue triangles). 
Inset: corresponding l/t avG vs. / showing a clear angular 
dependence of l/t ave for / — >> 0. 



ing to decoherence [18]. The corresponding decay time 
r c , which is the coherence time for m.2, was obtained 
by fitting T(r) to a function of the form e~ T l Tc . We 
can compare r c with the average dwell-time t ave in mode 
ni2, obtained by analyzing the instantaneous STO fre- 
quency [15] . Figure. [3] shows the central experimental 
result: r c and t ave as functions of current at 140°, 196°, 
and 220°. Both r c and t ave depend on field angle and cur- 
rent. At 196°, t ave exceeds r c by over an order of magni- 
tude for currents in a range about the threshold current. 
Here, mode- hopping is insignificant and does not limit 
coherence, the STO has a well-defined single mode, and 
we show below that NSTO theory applies here. For 140° 
and 220° , r c and t ave are approximately equal (within ex- 
perimental uncertainty in t ave [15 ]). Here, mode- hopping 
limits coherence time; more extensive analysis shows that 
for ip < 165° and cp > 205°, mode- hopping is the domi- 
nant decoherence process. 

We now discuss the experimental results. The mag- 
netization dynamics can be described by the Landau- 
Lifshitz-Gilbert (LLG) equation, which we write as 



dm A 

— = -mx[H 



-amx < m x 



Hn-^M 



- m x 



hd}, 



where m(r) is the local magnetization direction of the FL 
and M that of the RL; Ho is the total static effective field, 
including the out-of-plane spin-torque, hd is the dynamic 
demagnetizing field arising from the oscillating magne- 
tization density in FL, aj is the in-plane effective field 
due to spin-torque, and a the dimensionless damping. 
First, micromagnetic simulations confirm that the asym- 
metry in mode frequencies about ip = 180° [Fig. [lib)] can 
be caused by small ellipticity in the structure, with the 
exchange bias slightly misaligned with the major axis, 
as the resonance modes are sensitive to the equilibrium 
magnetization details. The magnetoresistance, however, 
remains nearly symmetric about cp = 180°, as it mea- 
sures the average angle between RL and FL magnetiza- 
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Figure 4. (color online). Wigner transforms of micromagnetic 
time-traces of average FL magnetization at T — K (left) 
and T — 300 K (middle). The right represents the Fourier 
transform (FFT) of the 100 ns-long time-trace at 300 K. The 
top row is for 140°; the bottom row for 180°. Bias current 
was 8 mA. 



tion directions. Now consider a system with the equi- 
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librium FL magnetization along Ho, which is at an an- 
gle (j) to the applied field. The linearized (about static 
equilibrium) LLG equation can be recast as an eigen- 
value equation £(r, r / ;cj)m(r / ) = 0. At zero damping 
and current, all torques in the LLG equation are con- 
servative, and the eigenvalues are real and undamped, 
but finite damping or current gives rise to nonconserva- 
tive torques, and eigenvalues become complex. We now 
outline how this leads to mode-coupling [15] . We restrict 
the discussion to two modes, and write the magnetization 
texture as m(r, t) = ci(t)mi(r)e~ %UJlt + C2(t)rh2(r)e~ %UJ2t , 
where rhi (r) is an eigenmode of the linearized LLG with 
eigenfrequency u^; the coefficients q carry the slow time 
dependence (on observational timescales r bs) of the mag- 
netization texture. We insert this into the full nonlinear 
LLG, project onto m^, and average over fast timescales, 
before expanding the nonlinear projections in q. We 
write the result as 



ci = -i [cji77 M |ci| 2 + ^ 2 r?i, 2 |c 2 | 2 ] ci -T G [1 + P m cji|ci| 2 + Pi, 2 o; 2 |c 2 | 2 ] c i + ^ [l -Qi,i^i|ci| 2 - Qi,2^2|c 2 | 2 ] Ci + i?i >2 C2 

C2 = ~i [uiV2,l\ci\ 2 +£J 2 77 2)2 |C 2 | 2 ] C 2 -T G [l+P 2) iCJi|ci| 2 +P 2 , 2 CJ 2 |C 2 | 2 ] C 2 +CT i [l -Q2,l^l|ci| 2 - ^2,2^2 |c 2 | 2 ] C 2 +i? 2 ,lCi, 



where rjij, Pij, and Qij are real; Rij is complex. The 
factors rjij are nonlinear frequency shifts, Pi j the nonlin- 
ear positive damping, Qij the nonlinear negative damp- 
ing, and <Jo the usual spin-torque coefficient defined in 
Ref. [3] These equations are generalizations of the equa- 
tion for the single- mode NSTO |3]. We note the lin- 
ear "backscattering" term i?i >2 (P 2 ,i) in these equations. 
This term does not appear in normal spin-wave expan- 
sions of Hamiltonians or equations of motion as it violates 
energy conservation on short timescales. Here, however, 
we consider conditions close to threshold: energy is ap- 
proximately conserved on long time-scales, and we en- 
force the constraint uj\ \ c\ | 2 +cj 2 1 c 2 1 2 = p, where p is a con- 
stant. On general grounds, we can set Qi, 2 = Q 2 ,i = Qo, 
Pi, 2 = P 2 ,i = Po, and i?i j2 = P 2 ,i = u + iv. We in- 
troduce new variables Qi^ 1 = y/uHCi and ijj = 2 — <pi 
and use the constraint QiQi + Q2Q2 — 0- The equa- 
tions for slow time evolution can then be recast as a two- 
dimensional Z 2 -invariant dynamical system, analogous to 
that describing SRLs p~9j[20]. Close to threshold, these 
equations are known [20] to have two stable solutions (ei- 
ther of the two modes) close to a homoclinic bifurcation. 
In the presence of a stochastic field originating from, e.g. 
contact with a thermal bath, competition between linear 
and nonlinear couplings can lead to mode-hopping [T9l - 
|2T] between the two stable solutions. Nonconservative 
torques lead to linear and nonlinear couplings between 
modes, as expressed by the terms in Qij, Pij, and Rij, 
i y^ j. As the bias current increases from zero, the local 
nonconservative field is m x [Hq — ^-M] ; this field also 
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dominates the damping (imaginary parts of eigenvalues) . 
For arbitrary field angles, a current exists — the thresh- 
old current-at which the energy dissipation rate equals 
the rate at which energy is pumped into the system by 
the spin torque. At this current, energy conservation 
Qi+Q 2 = is strict, but the couplings between modes do 
not vanish, and the system exhibits mode-hopping: the 
system's total energy is conserved, but energy is trans- 
ferred back and forth between the modes. However, at 
180°, Ho and M are collinear, and at threshold energy 
dissipation equals energy pumping and local dissipative 
torques vanish. Ignoring the nonlocal dissipative torques, 
this means that all terms multiplying Tq, cr 7, as well as 
Rij, vanish. Consequently, mode-hopping is minimized. 
The system's total energy is conserved, and each mode's 
energy is also individually conserved. 

Figure. [4] shows Wigner transforms of magnetization 
time-traces obtained from micromagnetic modeling for a 
system similar to the experiment [15 . At a finite tem- 
perature of T = 300 K and 1 = 8 mA, mode-hopping 
occurs at 140°, while at 180° precession is coherent at 
the m 2 mode — this current is just above the threshold 
for this field angle. Yet at T = K, no mode- hopping 
occurs. Apart from demonstrating that mode-hopping 
exists within the micromagnetic model, the modeling 
shows that (a) mode-hopping is induced by thermal fluc- 
tuations, and (b) the angle-dependence of mode-hopping 
agrees qualitatively with our analysis above. 

If the mode resident time t ave is long enough, one can 
analyze the system in terms of single-mode NSTO the- 



ory |3]. The experimental coherence time r c can then 
be compared with that from NSTO theory. Under the 
assumption that r c ex 1/A/, with Af the linewidth, we 
plot 1/(7tt c ) against bias current and fit 1/(ttt c ) to the 
NSTO expression [3 for subthreshold linewidth, A/ = 
Tg(1 — I /Ith), where Tq is the natural FMR linewidth, 
and / t h the threshold current [Fig.[5ja)]. For 196° we ob- 
tain r G « 300 MHz and I th = 6.4 mA. This T G is compa- 
rable to previous reports [22] and agrees well with an esti- 
mate using these material parameters for FL: saturation 
magnetization Mo = 1000 emu/cm 3 , Gilbert damping 
parameter cvq — 0.01. Above 7 t h = 6.4 mA, l/(7rr c ) in- 
creases with current, qualitative agreeing with the NSTO 
theory prediction. Thus the linewidth from NSTO the- 
ory qualitative describes 1/(ttt c ) at 196° well [15], con- 
sistent with £ ave being large enough that the STO is well 
described as a single-mode oscillator with decoherence 
caused by single- mode thermal fluctuations. However for 
140° and 220° the coherence is now limited by mode- 
hopping, and we show that the strong increase of 1/(ttt c ) 
above 8 mA is primarily due to increased mode-hopping, 
not due to increased nonlinearity in NSTO theory. 

Figure p\ (c) shows the measured power restoration 
rate T p [23 . According to NSTO theory, T p vanishes 
as rc(///th — 1) near the threshold. Again, for 196° 
the agreement between the measured power restoration 
rate and the NSTO theory prediction is reasonably good. 
But for 140° and 220° the power restoration rates, like 
the linewidths, have much larger minimum values, consis- 
tent with a limiting decoherence process other than ther- 
mal fluctuations about ni2. We note that even for 196°, 
T p does not vanish at / = i" t h? indicating either nonzero 
mode-hopping or the presence of additional sources of 
noise. Further inconsistencies between NSTO theory and 
experimental data is provided by the power distribution 
functions shown on the right of Fig. [5] The current values 
are chosen such that I // t h = 1-25 in all cases. According 
to NSTO theory, the power distribution function V(p) 
has the form of a Gaussian [24] exp[— (p — po) 2 /2Ap 2 ] for 
I > Ah, where po denotes stationary power and Ap the 
power fluctuations. This agrees well with the measured 
power distribution function near ip = 196°, shown by the 
solid line in Fig.[5|d), with the relative power fluctuations 
Ap/po in the 0.2-0.4 range, similar to that obtained in 
Ref. [25 j. In contrast, the power distributions are expo- 
nential at 140° and 220° , consistent with below-threshold 
conditions. We conclude that the threshold currents for 
these angles extracted from the linewidth measurements 
using NSTO theory are too small. Consequently, the in- 
crease in measured linewidth l/(7rr c ) above 8 mA at 140° 
and 220° cannot be attributed to an increase in nonlin- 
earity, unlike Ref. 22. This is further supported by the 
fact that the experimental nonlinear frequency shift did 
not show any strong angular dependence [T5] . 

One can qualitatively understand how moderate mode- 
hopping changes NSTO theory. First, it leads to effec- 
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Figure 5. (color online), (a) 1/(ttt c ) vs. current at 140° 
(black squares), 196° (red circles), and 220° (blue triangles). 
Solid lines are fits to the subthreshold region; Intercepts with 
the horizontal axis indicate threshold currents, (c) Power 
restoration rate T p vs. current at 140° (black squares), 196° 
(red circles), and 220° (blue triangles). Solid lines are visual 
guides. Right: Normalized power distribution function and 
fits at (b) 140°, (d) 196°, and (e) 220°. 



tively increased power dissipation from mode ni2, since 
hopping decreases the power in this mode. Second, for 
short excursions to other modes, the oscillator phase is 
random as it returns to ni2 , and is not correlated with the 
phase before the excursion. This means mode-hopping 
adds extra phase noise to oscillators. Within this sim- 
ple picture, we write the expression for the subthreshold 
linewidth as Af « Tq — cfqI + ^^; a is a dimension- 
less constant of order unity. Since a/t ave depends on the 
angle for / — » (see Fig. [3j inset), this explains the ob- 
tained apparent angular dependence of the zero-current 
linewidth Tq [Fig. [5ja)]. Furthermore, the extrapolated 
threshold current value shifts to a larger value, the shift 
being larger for weaker current-dependence of t ave , con- 
sistent with a larger threshold current for 140° and 220° 
than that extrapolated using NSTO theory. 

In summary, we have observed and analyzed multi- 
mode excitations and mode-hopping in an MgO-based 
spin-torque oscillator (STO). Mode-hopping occurred at 
all angles and currents, contrary to two-mode theory [5 
and expectations [3]. Insofar as the basic physics of 
mode-hopping appears to be intrinsic to the LLG equa- 
tions governing the dynamics of the oscillator nano-sized 
systems with single-domain equilibrium magnetization, 
efforts to reduce decoherence may need to focus on bet- 
ter understanding the energy barrier separating stable 
modes. 
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